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BAND STRUCTURE OF W AND Mo BY THE EMPIRICAL 
PSEUDOPOTENTIAL METHOD 
C. Guha Srldhar* and Ellis E. Whiting 
Ames Research Center 


INTRODUCTION 


The purpose of this paper is to present a calculation of the band struc- 
ture of W and Mo by the empirical pseudopotential method (EPM) . Both of 
these elements are in group VI of the periodic table, are transition metals, 
have a body centered cubic crystal structure and have interesting catalytic 
and superconducting properties. The availability of reliable band structures 
of these transition metals would be very helpful in understanding, theoreti- 
cally, many of their physical and chemical properties and to help predict the 
properties of suitably chosen alloys. 

The calculations of the band structure of W and Mo made previous to the 
present study are summarized in the next section. Although these calcula- 
tions have been very valuable, there is some question about the accuracy of 
the approximation used and it is always helpful to have other independent 
calculations for comparison. 

Recently, the optical reflectivity of both W and Mo was measured by 
Weaver, Lynch and Olson (refs. 1 and 2) over a wide frequency range. The 
availability of these data made it possible for us to undertake a new and 
independent calculation of band structure of W and Mo by using the empirical 
pseudopotential method. 

A complete description of the EPM is given in reference 3 and a few 
appropriate comments are presented in the Method section. However, as a 
means of introduction, in the EPM the crystal potential is approximated by an 
expression containing several adjustable parameters, called form factors. An 
approximate set of form factors are chosen to initiate an iterative procedure, 
A band structure is calculated using these form factors and then the band 
structure is used to calculate those properties that are to be used to indi- 
cate empirical adjustments to the form factors. After the form factors are 
adjusted, to provide better agreement with the chosen properties, the process 
is repeated until the calculated properties are believed to be adequate. 

In the present study, the form factors were adjusted by comparing the 
calculated reflectivity spectrum with the data (refs. 1 and 2) mentioned 
above and the location of the calculated Fermi surface (refs. 4 and 5), and 
results are also given in references 6 and 7 for W and reference 8 for Mo, 
In addition, the calculated charge density in the (110) plane was plotted 
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and used to indicate adjustments to the form factors to give a reasonable 
qualitative appearance to the charge distribution to obtain a stable crystal. 


PREVIOUS CALCULATIONS 


TV- The first band structure calculation for tungsten was done by Manning 
and Chodorow (ref. 9), using a cellular method. Later Loucks (ref. 10) (for 
the Cr group metals) and Matthiess (ref. 4) (for W) independently calculated 
the band structure and Fermi surface by using the nonrelativistic augmented 
plane wave (APW) method. More recently Petroff and Viswanathan (ref, 11) 
calculated the band structure and photoemission spectrum of tungsten by the 
APW method. Also, Lomer (ref. 5) has proposed a model Fermi surface for the 
chromium group metals based on Woods APW calculation on iron (ref. 12). 

Matthiess (ref. 13) has adjusted his band structure by introducing the 
spin orbit coupling to fit the Fermi surface data. Loucks (ref, 34) later 
modified his band structure calculation by introducing relativistic correc- 
tions in the APW calculation and then calculated the Fermi surface. More 
recently, band calculations by the Relativistic Augmented Plane Wave (RAPW) 
method and experimental volume and surface photoemission spectra have been 
published by Christensen and Feuerbacher (ref. 15). 

Mo- The band structure and Fermi surface of Mo were calculated by Loucks 
(ref. 10). Petroff and Viswanathan (ref. 11) have calculated the band struc- 
ture and density of states. Recently Iverson and Hodges (ref. 16) derived 
nonrelativistic energy bands by the renormalized atom method and determined 
the Fermi surface using an interpolated band structure. They also used the 
d e-Haas -- van- Alphen data for the neck of the electron "jack" and the electron 
"lens" to derive the spin orbit parameter. 

The density of states was calculated by Geguzin et al. (ref. 17) by using 
the Green's function method. Pickett and Allen (ref. 18) used the Slater 
Roster interpolation scheme on the APW band calculation of reference 11 to 
find the density of states. Koelling et al. (refs. 19 and 20) have calcu- 
lated the optical properties, the Fermi surface and the density of states. 


METHOD 


The diagonalization of large secular determinants in the pseudopotential 
calculations involves large amounts of computer time. Thus, an efficient 
iterative method for searching for an empirical pseudopotential which provides 
a near optimum agreement with experimental data is important. In the past, 
an inspection method (ref. 3) and nonlinear optimization schemes have been 
used. These schemes had great success in the calculation of band structures 
where single k-points in the Brillouin zone play a dominant role in deter- 
mining the optical properties. However, there are many elements which display 
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very strong "volume" effects where the main peaks in the reflectivity spec- 
trum arise not from "critical points" but from large regions throughout the 
Brillouin zone. Any attempt to label a given peak in the reflectivity spec- 
trum by assigning a single energy separation at a single k-point would, thus, 
be useless. To account for this volume effect, we have developed a new effi- 
cient iterative technique for finding the near optimum pseudopotential. The 
details of this method have been described in reference 22. The five pseudo- 
potential parameters for W and Mo are given in table 1. 


RESULTS 


Band Structure 

Figures 1(a) and 1(b) show the band structures of W and Mo. Since these 
band structures are very similar, they will be discussed together. There are 
several striking differences between the APW calculations (refs. 4 and 11) 
and the present calculation. 

(i) A: Near its end points close to T, the present Aj band lies above 

the A 2 band, but dips below this band at a point halfway between the F 
and H symmetry points. No dip is observed in the work of Matthiess (ref. 4) 
and Petroff and Viswanathan (ref. 11), in which the A^ band has a free 
electronlike shape and lies well above the A 2 and A 5 bands. 

(il) E: The N 4 state is low compared to the APW (refs. 4 and 11) 
results. As a consequence, for W the £4 band crosses the Ei band. 

In addition, for both W and Mo the present E 4 band is relatively parallel 
to the E 2 band over most of the plot. Thus, even though the dipole matrix 
elements for this band pair are small, the effect on the calculated reflec- 
tivity spectrum will be significant because of the large joint density of 
states associated with this parallel band pair. 

(iii) The present band ordering at the H symmetry point is somewhat • 
unusual. That is, Hj lies below H 15 , while one expects the Hi (7S) state 
to be higher in energy than Hi 5 ( 6 p). The expected ordering is found in 
references 4 and 11. In EPM calculations, the higher level (>5 eV) energy 
states are inherently less accurate than the lower energy states because of 
the truncation of the form factors at lower wave numbers. The assumed correct 
ordering might be achieved by using a nonlocal p-potential. However, we have 
avoided any such alteration of the pseudopotential because the main intent of 
the present work is to provide an accurate description of the bands lying 
within about 5 eV of the Fermi level, and also to keep the number of adjust- 
able parameters to a minimum. 

(iv) In the case of Mo, the Pi and P 3 states are reversed in compari- 
son to W. Also, the separation Pi - P 3 is larger for Mo than for W. 

In the APW calculation (refs. 4 and 11), the Pi state lies above the P 3 
state. 
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(v) Table 2 shows the comparison of some band separations. For W, the 
s-d separation is larger than the previous calculation (ref. 4) and the 
occupied d-bandwidth is less than the previous calculation. The s-p and 
d-bandwidths seem to be in reasonable agreement with the previous results. 

(vi) In the present calculations, T 2 s' ** Fi for Mo has a value that 
is the average of the values found in the two previous calculations (refs. 11 
and 16). However, Hgci - Fj has a higher value compared to both earlier 
calculations. The value of the d-bandwidth and occupied d-bandwidth from the 
present calculations both have smaller values than the earlier results 
(refs. 11 and 16). The occupied d-bandwidth (Ep - H 12 ) is especially very 
much smaller than the previous results. On the other hand, the present 
occupied s-d bandwidth (Ep - Tj) is larger than the two earlier results. The 
origin of the differences mentioned above could be due to a stronger s-d 
hybridization effect in the EPM, and the truncation of the local pseudo- 
potential, which would cause the higher s-p bands to be too low. 


Density of States 

The calculated density of states, N(E) is given in figures 2(a) and (b) 
for W and Mo, respectively. The histograms are calculated with 0.1 eV 
resolution. The cutoff point is the Fermi energy and the estimated uncer- 
tainty of its position is about 0.05 eV. For W, the theoretical result is 
superimposed on the experiment of Zeisse (ref. 23). Table 3 lists the peak 
locations of the density of states for the present calculation and for 
previous calculations and experiments. 

For W, the agreement of the density of occupied energy states is good 
to within about 0.3 eV with both previous APW calculations (refs. 4 and 11). 

The density of states at the Fermi level N(Ep) is 0.3 states of one spin/eV 
atom and is in good agreement with the APW calculations. However, McMillan 
(ref. 24) obtained a value of 0.15 states of one spin/eV-atom from experimental 
values of the superconducting transition temperature and the electronic speci- 
fic heat by using a strong coupling formulation. The agreement with McMillan’s 
result is poor. 

For Mo, the two experimental photoemission measurements show DOS peaks 
below the Fermi level of -0.5, -1.6, -3.9 eV from reference 25 and of -1.5 
and -3.7 eV from reference 26. Geguzln et al. (ref. 17) calculated DOS by 
using the Green’s function method and found peaks at -1.5, -2.7, -3.0, -3.8 
and -4.2 eV, Petroff and Viswanathan (ref. 11) calculated the same peaks by 
the APW method. In the present calculation, peaks are found at -0.5, -1.5, 
-2.5, -3.15 and -4.1 eV. It is clear that none of the calculations or the 
measurements gives completely consistent results. Thus, it is difficult to 
conclude whether our result is more or less accurate than the others. How- 
ever, the results of the density of states at the Fermi level agree very 
well. McMillan's (ref. 24) strong coupling formulation gives 0.28 states of 
one spin/eV-atom, as does Pickett and Allen (ref. 18), Petroff and Viswanathan 
(ref. 11), and the present study, whereas Geguzin et al. (ref. 17) finds a 
value of 0.31 spin/eV-atom. 
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Reflectivity 

In this subsection and the next two, we show the final Iterated value?,) 
of the calculated properties that were used to adjust the form factors. The 
good agreement between the final calculated and experimental results indicates 
that the resultant band structure discussed above is reasonably accurate. 

The calculated and experimental spectra are shown in figures 3 and 4. 

There are two theoretical curves on each figure, one for the nondirect tran- 
sition method. The iteration scheme used only the spectrum given by the 
direct transition model to adjust the form factors. After convergence, the 
curve for the nondirect transition method was calculated using the same form 
factors. The direct transition model predicts a spectrum which has a reflec- 
tivity approximately 10% higher than the data. Although this agreement is 
considered very good, the difference may be in part due to the relative 
inaccuracies involved in using the pseudo wavefunctlon tc calculate the 
dipole transition matrix elements. The nondirect model assumes the noncibn- 
servatlon of crystal momentum k. Both curves are shown because there has 
been significant discussion of both in the literature (refs. 27-25). 

Clearly the calculated reflectivity by the direct transition method is 
in quite good agreement for W up to 4.5 eV and for Mo up to 3.5 eV. 

Further work to improve the calculation at higher energy levels could be 
done, but the present results appear adequate for our purposes. 


Fermi Surface 

The central cross sections for the Fermi surfaces of W and Mo are 
shown iyi figures 5 and 6. The calculated Fermi energy contours are shown 
as the solid curves. 

For W, the experimental data (broken contours) are taken from Walsh and 
Grimes (ref. 6) and from Sparlin (ref. 7). The data of Walsh and Grimes 
give the linear dimensions of that portion of the electron "jack" and the 
hole octahedron in the (110) plane and the data of Sparlin give the hole 
ellipsoids at symmetry point N. Sparlin' s data were taken from his de-Haas- 
van-Alphen results. The agreement between the measured and calculated Fermi 
surface dimensions in the (110) plane is fairly good. The fact that the elec- 
tron "jack" and the hole "octahedron" fail to touch along A in experiment 
is readily explained in terms of the spin orbit coupling. This separation in 
W is estimated to be 5 % of TH distance. Matthiess and Watson (ref. 13) have 
estimated the 5d spin orbit coupling parameter in metallic W to be approxi- 
mately 0.4 eV. Spin orbit coupling also causes the "lenses" Inside the "necks" 
to reduce in size and perhaps disappear entirely. 

The dimension of the hole ellipsoids at N is determined by the energy 
difference between E(Ni') and the Fermi energy. Since Nj' is primarily a 
6-p state, this energy difference depends on the energy separation between 
the s-p and d bands. 
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The general features of the calculated Fermi surface of Mo are shown in 
figure 6 to be very similar to that obtained by Boiko et al. (ref. 8). The 
Fermi surface result agrees very well with Lomer's calculation (ref. 5) which 
is not shown in the figure. 

The size of the Fermi surface obtained experimentally by Boiko et al. 
(ref. 8), where in experiment the relativistic effect is naturally included, 
was smaller for W than for Mo. They also observed the vanishing of the 
electron lens in W. As far aB the present calculation is concerned, the 
Fermi surface dimensions of W are not significantly less than those for 
Mo. Precise values of the areas of Fermi surface have not been calculated 
because of the relative crudeness of the present interpolation scheme and 
because the calculation would Involve the comparatively inaccurate method of 
adding up the grid squares. 


Charge Distribution 

The charge distribution rras calculated on a mesh using 70 points in 1/48 
the Brillouln zone. The charge density for each valence band can be 
written as 

P n ® */d 3 k 2 |* nE (p E «| 2 

1 r 1 

where is the wavefunction for the nth band at k and {P^} is the 

symmetry operator needed to generate selected wave functions. The charge 
distribution of a completely filled band is normalized \ to two electrons/band 
in the primitive cells. The total valence band charge density is the sum 
over p n of the completely and partially occupied bands. 

Figures 7 through 10 show the individual valence bands and the sum of the 
valance bands charge densities. These charge densities are plotted on the 
(110) plane of the crystal. 

In each case, figure (a) is for W and figure (b) is for Mo. Since 
there are no experimental charge distribution results available for transi- 
tion metals, we will discuss the results of the two metals comparatively. 

In figure 7(a) for the first valence band of W, we observe some nonspher- 
ical features and lobes, whereas in figure 7(b) for Mo a nearly spherical 
symmetry of the charge distribution is seen at the atoms. This is because 
in Mo, the s-d separation (T i - H 12 ) in the lowest lying state is 3.3 eV, 
whereas in W it is only 1.9 eV. Since in W the d-like state Hj2 lies 
quite close to the s-like state, Tj, there is some hybridization which produces 
the lobes on the sides of the central atom. 

In the second valence band of W (fig. 8(a)) we do not observe any long 
lobes on the two sides of the central atom, as observed in Mo (fig. 8(b)), 
but instead we see small lobes along the body diagonal that indicates the 
contribution of d-llke states. The two long lobes on the sides of the central 
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contribution of d~like states. The two long lobes on the sides of the central 
atom arise from the P^ state of xy symmetry. 

The third valence band (figs. 9(a) and (b)) shows fairly strong maxima 
along the diagonal* which indicate major contributions from the d-states. 

The 4th and 5th valence bands are mostly empty below the Fermi level, and 
their charge distributions are not shown. 

The total charge distribution of the occupied states is shown in fig- 
ures 10(a) and 10(b). The maximum of the contours associated with each atom 
occurs near l/4th of the interatomic distance. These are made up primarily 
from the lobes of the d-states. In the case of W, the lobes are consider- 
ably smaller than those in Mo, although the magnitude of the contours is 
comparable. This difference may be due to the fact that (P^ - Hj 2 ) is only 
1.0 eV in Mo and 2.25 eV in W. 


SUMMARY 


> \ 

An effective new scheme to determine the empirical pseudopotential has 
been used to calculate the band structures of W and Mo. The band structures 
for the two metals are similar, as expected, and both differ In several sig- 
nificant ways from previous APW calculations. 

The present occupied density of states in W agrees very well with the 
previous APW calculations and experiment. The agreement in the case of Mo 
is fair with other experiments and theory. However, it is difficult to 
resolve the density of state peaks from photoemission measurements. Hence, 
the poor agreement, in the case of Mo with experiment, is not considered 
to be serious. The density of states at the Fermi surface shows excellent 
agreement with previous APW (refs. 4 and 11) calculations in the case of W. 
However, a reason for the poor agreement with McMillan's result (ref. 24) 
is unknown. In the case of Mo, this result shows good agreement with all 
the previous calculations (refs, 11, 17, 18, and 24). 

The reflectivity of W and Mo has been calculated from 0-4 .5 eV . The 
results look very good except for the region above 3.5 eV for Mo. The EPM 
iteration used the reflectivity spectrum calculated by the direct transition 
method to adjust the form factors. After the form factors were found in this 
way, a second reflectivity spectrum was also calculated using the nondirect 
transition method. The results from the two methods show some differences. 

The cross sections of the Fermi surfaces for W and Mo have been mapped 
out. They show close agreement with an available data. W shows good agree- 
ment with Matthiess result and Mo shows good agreement with the original 
Lomers (ref. 5) model Fermi surface. The charge density was calculated for 
the (110) plane and show a qualitatively correct charge distribution. We 
believe that the present band calculations of W and Mo are sufficiently 
accurate to be used to calculate phonon coupling constants and bonding proper- 
ties in alloys of these transition metals. 
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TABLE 1.- EPM PARAMETERS 3 



Present results, Ry 

Parameter 

Tungsten 

Molybdenum 

V (2) 

-0.050 

-0.082 

V(4) 

.035 

.101 

V(6) 

.015 

.048 

V (8) 

.109 

.069 

Depth of d-well 
Width of d-well 

-3.96 

-3.90 
1.09 A 


3 The argument of the form factors is 

given in units of (2ir/a) 2 where a = 3,16 A 

for W and a = 3.15 A for Mo. The value 

of the free-electron parameter is 

y = 0.22xi0 -1 * 9 g-cm 3 * * * * * for W and 

0. 320*10~ 49 for g-cm 3 for Mo, 
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TABLE 2.- INTERESTING BAND DIFFERENCES 
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TABLE 3.- COMPARISON tV DENSITY OF STATES 


o 








£ «. £ 
I I I ui Z 



A 8 'AOU3N3 


12 


POSITION IN BRILLOUIN ZONE 






S2 -fi 5! 

ill i 



A a 'ADH3N3 


13 


POSITION IN IRREDUCIBLE BRILLOUIN ZO*'E 




iiNn Aavunauv oin soa 


14 


Figure 2.- Density of states histogram. 
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Figure 2.- Concluded. 
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Figure 4.- Reflectivity spectrum of Mo. 
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Figure 5.- Central portion of Fermi surface of 
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Figure 7.- Concluded. 
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SUM OF VALENCE BANDS 



Figure 10.- Charge distribution of occupied states (sum of valence bands) 
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Figure 10.- Concluded. 


